Formation of a galaxy with a central black hole in the Lemaitre-Tolman model 



Andrzej Krasinski 
A''. Copernicus Astronomical Center, 
Polish Academy of Sciences, 
Bartycka 18, 00 716 Warszawa, Polan^ 



o 
o 

(N 

a 

(N 
> 



a\ ■ 
o : 

m . 
o 



Charles Hellaby 

Department of Mathematics and Applied Mathematics, 
University of Cape Town, Rondebosch 7701, South Afric^ 
(Dated: 25/9/2003) 

We construct two models of the formation a galaxy with a central black hole, starting from a 
small initial fluctuation at recombination. This is an application of previously developed methods 
to find a Lemaitre-Tolman model that evolves from a given initial density or velocity profile to a 
given final density profile. We show that the black hole itself could be either a collapsed object, or a 
non-vacuum generalisation of a full Schwarzschild-Kruskal-Szekeres wormhole. Particular attention 
is paid to the black hole's apparent and event horizons. 
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I. AIM AND MOTIVATION 

It has become generally accepted that most large 
galaxies contain central black holes (e.g. This 
is due to mounting evidence for very high luminosities 
within very small radii at the centre of our galaxy and 
many others, as well as high orbital velocities of stars 
very close to the centre, and is bolstered by observations 
of radio and optical jets, as well as Seyfert galaxies and 
quasars at large redshifts. Since the 'very small' radii ac- 
cessible by current observations are still well outside the 
horizons of the putative black holes, their exact nature, 
and even their existence, is open to debate. However, 
since the mean density inside the horizon of a black hole 
is ps = 3cV327rG3M2 = 1.845 x 10^'^ /{Ms/Mq)^ g/cm^, 
it is difficult to see how black hole collapse can be avoided 
above lO'^ Mq. 

The Lemaitre-Tolman (L-T) model describes the be- 
haviour of a spherically symmetric dust distribution and 
has been a very fruitful source of models of inhomoge- 
nous cosmology, smaller scale structure formation, and 
even black holes and naked singularities. 

In paper I we considered the problem of finding 
a spherically symmetric model that evolved from a given 
initial density profile to a given final density profile. We 
showed that this can always be done with an L-T model, 
and we developed an alogorithm to find the arbitrary 
fimctions of such an L-T model from the given profiles. 
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A numerical example produced an Abell cluster from a 
density fluctuation at recombination. 

In paper II [i^, we generalised to finding L-T models 
that evolve from a given velocity profile to a given density 
profile, the converse, and also between two velocity pro- 
files. Several numerical examples, including the evolution 
of a void, demostrated the usefulness of the method. 

We now utilise these methods to consider the formation 
of a galaxy with a central black hole, a task for which the 
L-T model is particularly well suited. Although spiral 
galaxies are not exactly spherically symmetric, both the 
core and the halo — together containing more mass than 
the disk — are quite close to it, so the L-T model is not 
a bad first approximation. 

The present state of the galaxy is defined by a mass 
distribution that consists of two parts: 

1. The part outside the apparent horizon at t2 ~ for 
which we use an approximation to the observationally 
determined density profile of the M87 galaxy. This part 
extends inward to a sphere of mass Mbh, where Mbh is 
the observationally determined mass of the black hole in 
the M87 galaxy. 

2. The part inside the apparent horizon at t2. Since, 
for fundamental reasons, no observational data at all ex- 
ist for this region apart from the value of Mbh, we were 
free to choose any geometry. We chose two examples: 

2a. A simple subcase of the L-T model, discussed in 
sec. Ill F as an illustrative example of properties of hori- 
zons. In this model, the black hole does not exist initially 
and is formed in the course of evolution. 

2b. A pre-existing wormhole, also chosen arbitrarily 
for simplicity of the calculations. 

The boundary between the "inside" and "outside" at 
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times other than t2 goes along a comoving mass sheh, 
so that at t < t2 the apparent horizon resides in the 
inside part. 

For the initial state, at ti = (the last scattering of 
the cosmic microwave background (CMB) radiation), no 
usable observational data are available, either, but hope- 
fully only temporarily. The expected angular size on the 
CMB sky of a perturbation that will develop into a sin- 
gle galaxy (0.004°) is much smaller than the current best 
resolution (0.2°). Therefore we tried an exactly homoge- 
neous initial density and a homogeneous initial velocity. 
The former turned out to lead to an unacceptable con- 
figuration at a collapsing hyperbolic model with no 
Big Bang in the past. Consequently, we settled on the 
homogeneous initial velocity, which then implied the am- 
plitude below 10^^ for the initial density perturbation. 

These two states, at ti and t2, uniquely define the L- 
T model that evolves between them, as shown in Paper 
II. The 3-d surface graphs of density as a function of 
mass and time show that the evolution proceeds without 
shell crossings, and so the model is acceptable, at least 
qualitatively. 



ii. basic properties of the 
lemaItre-tolman model 



The Lemaitre-Tolman (L-T) model [22,123 is a spher- 
ically symmetric nonstatic solution of the Einstein equa- 
tions with a dust source. See for an extensive list of 
properties and other work on this model. Its metric is: 



{R'f 



1 + 2E{r) 



- dr 2 - i?2 (t , r ) (di?2 -f sin^ Mip"^), 

(2.1) 



where E{r) is an arbitrary function of integration, R' 
dRjdr, and R{t,r) obeys 



(2.2) 



where 7? = dR/dt and A is the cosmological constant. 
Eq. I|2.2|l is a first integral of the Einstein equations, and 
M — M{r) is another arbitrary function of integration. 
The mass-density is: 



2M' 
R^R' 



where k 



SttG 



(2.3) 



In the following, we will assume A = 0. Then eq. (|2.2|) 
can be solved explicitly, and the solutions are: when E < 



(elliptic evolution): 
i?(t, r) 



pfgy(l-cosry). 



(2.4a) 



77-sin77 = ^ "^M ' (t-tBir)), 



M -vw;, ^ (2.4b) 

where 77 is a parameter; when E — (parabolic evolu- 
tion): 



R{t,r) 



-M{t--tB{r)f 



1/3 



and when E > Q (hyperbolic evolution) 



R{t,r) =i|(coshr;-l), 
sinhr;-77 = ■^^^(t - ^^(r)) 



(2.5) 



(2.6a) 
(2.6b) 



where tsir) is one more arbitrary integration function 
(the bang time). Note that all the formulae given so far 
are covariant under arbitrary coordinate transformations 
f = g{r), and so r can be chosen at will. This means 
one of the three functions E{r), M{r) and tB{r) can be 
fixed at our convenience by the appropriate choice of g. 
We can define a scale radius and a scale time for each 
worldline with 



P(r) 
Tir) 



2M 

W\ 

27rM 
|2£:|3/2 



(2.7) 
(2.8) 



and it is evident from (|2.4|) that, for the elliptic case, 
these are the maximum R and the lifetime for each r 
value. The crunch time is then 



tc{r)^tB{r)+T{r). 



(2.9) 



Writing eq. I)2.4b|l at r/ = 2tt, where t — tc, i.e. at 
the Big Crunch, and then dividing the two equations we 
obtain 



77 - sin 77 = 27r(t - tB)/{tc ~ ts), 



(2.10) 



so larger 77 means only that the dust particle has com- 
pleted a larger fraction of its lifetime between the Bang 
and the Crunch. 

The parametric solutions 1)2. 4|l and 1)2.6(1 can also be 
written 
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t^tB 



t = tB 




arccos | 1 H — j ^ \ 



TT + arccos 



-ER / ER^ 
M [ ^ IT , 



2ER\ „ -ER / ER 




(2.11) 
(2.12) 



r 



for the expanding and collapsing elliptic cases, and 



t^ts 



M 



(2£;)3/2 



I2ER ( 2ER\ 
'2 



M 



M ) 



arcosh 1 



^ER' 



(2.13) 



for the hyperbolic case (expanding). 

Apart from extended parabolic regions, there are also 
parabolic boundaries between elliptic and hyperbolic re- 
gions, where E ^ 0, but E' ^ 0. The limiting forms 
of equations H2.4|l and H2.6|l are found by requiring well 
behaved time evolution and setting 



rj = tj^/e 



(2.14) 



so that fj is finite if {t — Ib) is. 

The Friedmann models are contained in the L-T class 
as the limit: 



ts — const. 



|£:|3/2/Af 



const, 



(2.15) 



and one of the standard radial coordinates for the Fried- 
mann model results if the coordinates in 12. 4() - (|2.6|l are 
additionally chosen so that: 



M = Mnr^ 



E = EnV" 



(2.16) 



with Mq and Eq being constants. 

In constructing our galaxy model, it will be convenient 
to use M{r) as the radial coordinate (i.e. f = M(r)) 
— because in most sections we shall not need to pass 
through any "necks" or "bellies". Thus, M(r) will be a 
strictly growing function in the whole region under con- 
sideration. In some of the sections we shall consider a 
black hole with a "neck" or "wormhole" , but even there, 
because of spherical symmetry, we will consider only one 
side of the wormhole, where M(r) is also increasing. 

Then with R ^ R(t,M): 



up - 
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m OR 

dM 



dM 



(2.17) 



In the present paper we will apply the L-T model to 
a problem related to that considered in papers I & II: 



Connecting, by an L-T evolution, an initial state of the 
Universe, defined by a mass-density or velocity distribu- 
tion, to a final state defined by a density distribution that 
contains a black hole. 



A. Origin conditions 

An origin, or centre of spherical symmetry, occurs at 
r = rc ii R{t,rc) = for all t. The conditions for a regular 
centre have been derived in |25l | from the requirements 
that, away from the bang and crunch, and in the limit 
r — > rc'. 

u rj in H2.4|l and (|2.6|l must be finite if (t — t^) is 
finite, 

• the density (|2.3|) and the Kretschmann scalar are 
non-divergent, and the density is not zero, 

• on a constant time slice dp/ dR = 0. 

In the equations below, the symbol Od{M) will denote a 
function that has the property ImiM^o {Od{M)/M'^) = 
0. The resulting conditions for the neighbourhood of Vc 
are 

R = l3{t)M^/^ + Oi/3(M) along constant t, (2.18a) 



^ = 7M2/3 + 02/3(M), 
tB = tOc{M) , c > 1/3, 
up = 6/f5^ + Oo{M), 
M{rc) - 0. 

We also need r < to avoid shell crossings. 



(2.18b) 
(2.18c) 
(2.18d) 
(2.18e) 



B. Sliell crossings, maxima and minima 

Shell crossings, where a constant r shell collides with 
its neighbour, are loci of i?' = that are not regular 
maxima or minima of R. They create undesireable sin- 
gularities where the density diverges and changes sign. 
The conditions on the 3 arbitrary functions that ensure 
none be present anywhere in an L-T model, as well as 
those for regular maxima and minima in spatial sections, 
were given in [2^ . and will be used below. 
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III. APPARENT AND EVENT HORIZONS IN 
THE L-T MODEL 

We will be modelling a galactic black hole, so it will 
be useful to consider its horizons. Apparent and event 
horizons of L-T models were studied in |27l |. in which 
L-T models that generalise the Schwarzschild-Kruskal- 
Szekeres topology to non-vacuum were demonstrated. It 
was shown that, when there's matter present, the light 
rays get even less far through the wormhole than in the 
vacuum case. The diversity of possible topologies was 
discussed. We lay out further details of the apparent 
horizon here. 



A. Definitions and basic properties 



Let us write the evolution equation (|2.2|l with A = 



as 



2M 

I £ = +1 in the expanding phase, 
where < (3.1) 
I £ ^ —I in the collapsing phase. 



= ^dt^ + ^^dr\ 
1 + 2E 



The radial light rays must be geodesies by symmetry: 

(3.2) 

± -r ^J-^ 

and this may be written as 

jR' 



dt 

"~ dr 



„ VTT2E ' 

j = +1 for outgoing rays, 
j — —1 for incoming rays, 



where 



(3.3) 



whose solution we write a.s t = i,i(r, t„o)j or often just 
tnir) or tn. 



1. Apparent horizons 



Along a ray we have 

Rn = R(tn,r) 

(RnY ^ R^ 



R' 



1 R'. 



(3.4) 



(3.5) 



The apparent horizon (AH) is the hypersurface in space- 
time where the rays are momentarily at constant R: 



{Rn}' = 



2M 



2E = -£j\/l + 2E 



-1, 



and R = 2M. 



(3.6) 
(3.7) 



There are in fact two apparent horizons: 
The future AH: (AH+), where 

(outgoing rays) 
(in a collapsing phase), 

and the past AH: (AH^), where 




(3.8a) 




(3.8b) 



(incoming rays) 
(in an expanding phase). 

We find dt/dr along the AH by differentiating 1)3.7(1 : 

Rdt + R' dr ^2M' dr, (3.9) 



givmg 



, _ dt 



AH 



2M' - R' 
R 



2M' - R' 



2M 
R 



2E 



and, since R = 2M on the AH, 



dt 
dr 



AH 



£{2M' - R') 



(3.10) 



(3.11) 



In the vacuum case p = 0, which implies M' = 0, we have 
dt/dr I AH = dt/dr I „ since £j = -1. Note that M' = 
could be only local, so the AH would only be null in 
that region. In the Schwarzschild metric, where M' — 
everywhere, this is consistent with R = 2M being the 
locus of the event horizons; and in this case they coincide 
with the apparent horizons. 

Recall that in the Schwarzschild spacetime the future 
& past event horizons, EH+ & EH~, cross in the neck at 
the moment it is widest. (Call this event O.) This holds 
for L-T models too. For hyperbolic regions, with E > 
along each dust worldline, there is either only expansion 
or only collapse, i.e. only one AH (either AH"*" or AH^) 
can occur. The AHs can thus cross only in an elliptic 
-E < region. At the neck of a L-T wormhole, where 
2E = — 1, Af is a minimum, and is maximum, the 
moment of maximum expansion is 



i?2 







2M 
~R 



1 



Rmax{M^in) = 2Af. (3.12) 



At all other E values in an elliptic region — 1 < 2E < 0, 
we find 7? = R^ax = 2M/{~2E) > 2M. Thus 
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R = 2M has two solutions — one in the expanding phase 
& one in the collapsing phase. So the AH+ & AH~ meet 
at the neck maximum (event O). 

To establish whether AH is timelike, null or spacelike, 
we compare the slope of the AH+ with the outgoing light 
ray (or the AH~ with the incoming light ray), i.e. £j = 
-1: 



dr 



AH/ 



dt 

dr 



-ej 1 



2M' 



2M' 
(3.13) 



but below, we will actually calculate 



B 



M' 



AH± 

l + B 
1-B 



AH^ 



1-WI + 2E (3.14) 
dM 



where we have used (|3.11|) and written R' /M' — dR/dM 
and t'g /M' — dts/dAI, since M' > 0. Now since the con- 
ditions for no shell crossings require Af > where 
R' > and vice-versa, we have 



Bmax = 1, B = +00 AH+ outgoing null (when M' — , p — 0), 

1 > B > -1, +00 >B>0^ AH+ spacelike (for most M'), . . 

B = -l, 5 = AH+ incoming null (for large M'/R'), ^ ' 

-I > B > — oo, > -B > — 1 AH+ incoming timclikc (for very large M'/R'), 



so an outgoing timelike AH"*" is not possible. This means 
outgoing light rays that reach the AH"*" always fall inside 
AH+, except where M' = 0, in which case they move 
along it. The possibility that AH* is timelike {B < 0) 
holds if^ 



diAH 

dM 



1 



VT+2E 



(3.16) 



This means that the smaller E is, the steeper the locus 
of the AH must be to make it timelike. 

The argument is similar for light rays at AH^ , except 
that 'incoming' should swopped with 'outgoing'. If ingo- 
ing light rays reach the AH~ , they pass out of it or run 
along it. 

Since this is true for every point on AH+ & AH~ , radial 
light rays in dense wormholes are more trapped and go 
even less far than in vacuum. In particular, if^ p > 
where 2E = — 1, light rays starting at O fall inside AH+ 
to the future and inside AH~ to the past. 



2. Event horizons 

The event horizon is the very last ray to reach future 
null infinity (EH"^), or the very first one to come in from 
past null infinity (EH~). If we have vacuum {M' = 0) 
everywhere, then light rays travel along R = 2M, and the 
EHs coincide with the AHs. If there is matter, M' > 0, 
on any worldline, then the incoming light rays emerge 
from AH~ and outgoing light rays fall into AH+ at that 
r value, and so the EHs split off from the AHs (see ^7}). 



B. Locating the apparent horizons in elliptic 
regions 

1. AH~ during expansion 

We shall first consider the expansion phase of an el- 
liptic model, where < 77 < tt and i? < 0, so we have 
£ = +1 and only AH^ is present. Since R — 2M on an 
AH, we have from H2.4a|l : 



^ Note that eq. 13.161 can be equivalently written as follows: 

M' 



> 



V1 + 2-E 



■■M', 



where A4 is the sum of all the rest masses within 
the r = const sphere, equal to 47r y^—git, x)dx = 
4nJ^ [R^{t,x)R'{t,x)/V2E + 1] dx, while M 
4n jriR? {t, x)R' {t, x)dx is the active gravitational mass 
fsee l28l for a detailed discussion). 

Since a regular minimum or maximum _R' = requires all of M' , 
E' and t'g to be locally zero, p > ^ M"/R" > 0. 



C0S77AH = 1 + 4i?, 



(3.17) 



and thus, along a given worldline, the proper time of 
passing through the AH, counted from the Bang time 
tB, can be calculated from H2.11|l with R = 2M to be 



ts^M 



arccos (1 + 4£;) - 2^-2E{l + 2E) 



i-2E) 



3/2 



(3.18) 
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The function F = ^ah- ~ of the argument / = 2E, 
defined in H3.18|l has the foUowing properties 



Putting and (|2.18c|l in we find^, neglecting 

powers of M that are necessarily positive, 



Fi-l) 
dF 
df 



^ Mtt, F{0) = 4A//3, 
< for - 1< / < 0, 



(3.19) 



i.e. it is decreasing. These properties mean that the AH 
does not touch the Big Bang anywhere except at a centre, 
M = 0, even ii E = 0. 

Along all worldlines with E > —1/2, even parabolic 
worldlines E = Sz M ^ 0, the dust particles emerge 
from AH^ a finite time after the Big Bang (77 = 0) and a 
finite time before maximum expansion (77 = vr). In order 
to find the slope of AH~, we differentiate (|3.18|l with 
respect to M to obtain 



dt 



AH- 



dM 



1 



(-2£;)3/2 



-vr 

E 



2E - 



3M dE' 

(-2£;)5/2 dM_ 
Af (3 + 2E) 



arccos (1 + AE) 



dE 



2E^^/T+2E dM 



dtB 
dM' 



(3.20) 



In general, this is very difficult to analyse, but for special 
cases it will be possible. 

So although (|3.17|) shows wordlines with larger E exit 
AH~ at a later stage of evolution, this may not corre- 
spond to a later time t, or even to a longer time (t — ts) 
since the the bang. It is not at all necessary that i? is a 
monotonically decreasing function of r in an elliptic re- 
gion; in general it can increase and decrease again any 
number of times. 



dt 



AH- 



dM 



ctM" 



3' 



(3.23) 



and from H3.14|l with £ — +1 and H3.23|l we sec that. 



B 



-ctM" 



(3.24) 



The behaviour of the apparent horizon in the vicinity of 
the center is not unique, each of the cases listed in (|3.15|) 
can occur, depending on whether c > 1. This nonunique- 
ness of behaviour is connected with the shell-focussing 
singularities that appear in some L-T models. Various 
studies |2^-[33 have shown that outgoing light rays may 
emerge from the central point of the Big Crunch and even 
reach infinity. Moreover, even though, in those cases, this 
central point (where the Big Crunch first forms) appears 
to be a single point in comoving coordinates, it is in fact 
a finite segment of a null line in the Penrose diagram, 
see jl^ . At the Big Bang, we have the reverse — incom- 
ing light rays may reach the central point of the Bang 
singularity. Any radial light ray emitted from the cen- 
ter of symmetry that falls into the Big Crunch must first 
increase its R value, and then decrease, i.e. must cross 
the AH-I- in between. In consequence, the AH-I- cannot 
touch the center of BC earlier than the null singularity 
does. Similarly, the AH^ cannot touch the center of BB 
later than the null singularity does. 



2. AH near an origin 



On the AH in the neighbourhood of a regular center, 
we obtain from H3.17|l and (|2.18b|) to lowest order: 

1 - cosr/AH = -47M2/3 4- 02/3(M). (3.21) 

and consequently, ?7ah- ^ as r ^ Tc, i.e. AH~ touches 
the Big Bang set at the center. Since E{rc) — and 
-E < in the neighbourhood of the center, it follows that 
E'{rc) < 0, and, via H3.17|l . that the dust particles with 
larger r (smaller E) exit AH~ with larger values of rj. 

To show the behaviour of AH^ near rc is not unique, 
we take the following example for E{M) 

^ = M2/3(7 + 72Af'^) , 

(i>2/3, 7^0, 72 7^0. (3.22) 



3. AH in the parabolic limit 

A shell of parabolic worldlines occurs at the boundary 
between elliptic and hyperbolic regions, where E ^ 
but E' ^ and^ M > 0. From (|XT7|) . JSHJ, 



^ Since M^^^ is a natural measure of proper radius near an origin, 
the slope of AH— is zero: 

^^AH_ ~ ScrM'^-l/S + 4M2/3 + ... 0. 
dAfi/3 

Away from the centre, the sign of dt^jj- /dA/l/^ depends on 
whether or not c > 1 , as r < for no shell crossings ,2fS] . 
^ Elliptic and parabolic regions are not possible for A/ = 



7 



and lIXTIjl with £ 



-1, we see that 



l-AH 



tB 

T 

dM 



E^O 





4M 

oo 



E^O 



AM dE 

1 AM dE 
3 ^ 5 dM 



dtB 
dM 

_ dtB 

dM 



(3.25a) 

(3.25b) 
(3.25c) 

(3.25d) 

(3.25e) 



so AH" never touches the bang here, despite 77 being zero. 
The divergence of the worldhne hfetime T suggests that 
either the bang time or the crunch time recedes to infin- 
ity, as would be expected in a hyperboHc region. There is 
in fact a third possibihty, as there is no reason why both 
times should not diverge. Indeed, an asymptotically flat 
model is achieved by letting £' — > as r — > 00 [i^l, with 
both bang and crunch times diverging^. We also see the 
slope and the causal nature of AH" are uncertain here. 



4- AH^ during collapse 



0, tc - tAH+ 
d^AH+ 

dM 



AM/ 3, and 



E^O 



E^Q 



AM dE 
~dM 
AM dE 
~dM 



dtc 
dM' 
dtc 
dM' 



(3.28a) 
(3.28b) 



where tc is given by H2.9|l & H2.8|) . 



C. Apparent horizons in parabolic regions 

The corresponding results in an expanding E = 0, 
E' ~ L-T model follow from section HHB 31 (or directly 
from (jTsJ with R = 2M). In particular (|3.25dl) and 
(|3.25e|l become 



dt 



AH- 



dM 



4 dtB 

3 dM 

1 dtB 

3 dM 



(3.29a) 
(3.29b) 



Despite the no shell crossing condition d^s/dM < 0, 
AH~ may still exhibit all possible behaviours of (|3.15|) 
with "outgoing" and "incoming" interchanged. In a col- 
lapsing parabolic model, time reversed results apply. 



This can be obtained form the above by replacing (t — 
ts) with {tc — t), rj with {2tt — ry), flipping the signs 
of £ and j, and swopping "incoming" with "outgoing". 
However, keeping tB as our arbitrary function, we have 
TT < 77 < 27r. Eq. l|TT7|l still applies, but instead of (|XTH|) 
and H3.20|l . we now obtain 

{t - iB)AH+ = 

TT + arccos(-l - AE) + 2^ -2E{\ -h 2E) (3-26) 



M 

d^AH+ 

dM 



1 



(-2£;)3/2 



3M dE 



[-2Ef/^ {~2Ef/^dM 



[tt -f arccos (-1 - AE)] 



lVTT2E + ^^i^±3^^ + ^. 

E^ 2E^VTT2EdM dM 



(3.27) 



and H3.14|l applies with £ — —1. The special cases all 
follow in the same way. Near an origin, we again find 
nonunique behaviour, one of the possibilities this time 
being light rays escaping the crunch at rc before AH+ 
forms. At a regular extremum where E ~ —1/2, ?7ah+ ~* 
TT and AH+ crosses AH" . In the parabolic limit, 7?ah+ ^ 



D. Locating apparent horizons in hyperbolic 
regions 

Using the same methods as for expanding elliptic re- 
gions, we find the behaviour of the AH in expanding hy- 
perbolic regions is qualitatively the same. There is of 
course only one AH, no maximum expansion, and loci 
where E = —1/2 are not possible. But the results for 
origins and for the parabolic limit both carry over. Col- 
lapsing hyperbolic regions are essentially like collapsing 
elliptic regions. 



E. Locating the event horizon 

In general, locating the EH involves integrating (|3.3|l . 
which cannot be done analytically because R' is not a 
simple expression. In addition, we must find the last 
outgoing ray to escape the crunch, or the first incoming 
ray to avoid the bang. For a particular model, EH+ 
is located numerically using a compactified coordinate 
representation of the spacetime in section IIII F 31 



F. An illustration — a simple recollapsing model 



^ Another example of this kind will be discussed in sec. 1111 1^1 



We shall illustrate several properties of the L-T model 
and of its apparent horizons on a simple example. In the 
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present subsection, we will use this model with unrealistic 
parameter values, chosen in such a way that all the figures 
are easily readable. Later, in sec. IIVI we will use the 
same model for modelling a galaxy with a black hole at 
the center, with parameters chosen to fit observational 
data. 



The evolution of our model can be visualised more 
graphically in a 3-d diagram (Fig. |2Jl, which shows the 
value of the areal radius R at each time t and each value 
of AI. Fig. O shows the "topographic map" of the sur- 
face from Fig. [3 It contains contours^ of constant R (the 



1. Definition 



We take an i? < L-T model with a regular centre, 
whose Big Bang function tsiM) is 



tB{M)^~hM^ + tBQ, 

and whose Big Crunch function is 



(3.30) 



(3.31) 



where the parameter Tg is the lifetime of the central 
worldline where M — Q. The numerical values of the pa- 
rameters used in the figures will be a = 2 • 10**, b = 200, 
^so = 5, To = 0.05. Their values were chosen so as to 
make the figures readable and illustrative, at this point 
they are unrelated to any astrophysical quantities. Since 
t = tc a,t 1] — 27r, we find from cq. (|2.4|l : 



E{M) 



1 f 2t:M 
2 



tc — ts 

2\ 1/3 



2/3 



M2/3 



(aM3 + bAP + To 



,2/3 



(3.32) 



Note that as M ^ oo, we have ts — oo, tc +oo 
and E ^ 0. Hence, the space contains infinite mass and 
has infinite volume. Unlike in the Friedmann models, 
positive space curvature does not imp ly fi nite volume; 
this has been known since long ago |26ll38| . 



2. Description 

The main features of this model are shown in Fig. ^ 
Note that the AH"*" first appears not at the center, but 
at a finite distance from the center, where the function 
tAH+{M) has its minimum, and at a time ths < tc{0)- 

At all times after the crunch first forms, t > tc{Q), 
the mass Ms already swallowed up by the singularity is 
necessarily smaller than the mass Mbh that disappeared 
into the AH+. So for an object in which the black hole 
already exists, the time t — now must be taken after t^^g 
and in Fig. ^it is also greater than ic(0). The mass Ms 
cannot even be estimated by astronomical methods. 

The situation is similar, but reversed in time, for the 
Big Bang singularity and the AH~. 



crunch/ 



AH+ 



ths" 



center 
(M = 0) 



Ms; 



Mbh 



bang--,^ 



FIG. 1: Evolution leading to a black hole in the E < Q L-T 
model of eqs. 13.3011 - 13.3111 . The final state is defined at 
the instant t — t2 — now. Worldlines of dust particles are 
vertical straight lines, each has a constant mass-coordinate. 
Intersections of the line t — t2 with the lines representing the 
Big Crunch and the future apparent horizon determine the 
masses Ms and Mbh, respectively. 



^ Fig. 121 can be used to explain some properties of shell crossings. 
Since the worldlines of the dust source are constant M lines (i.e. 
vertical straight lines), they can never intersect in a t-M diagram, 
even at shell crossings. A shell crossing would show in this picture 
as a point where a constant- ij-contour has a horizontal tangent. 
(No such points exist in this case, because the functions of the L- 
T model had been chosen appropriately.) The figure graphically 
explains why, for avoiding shell crossings away from a neck or 
belly, it is necessary that the Big Bang is a decreasing function 
of M, and the Big Crunch is an increasing function of M. The 
center of symmetry, the Big Crunch and the Big Bang together 
form the R = contour. Contours of small constant R must 
have a similar shape. Hence, if either of the two conditions were 
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thinner curves) inscribed into Fig. ^ (The other curves 
are outgoing radial null geodesies, and we shall discuss 
them further below.) 

Rembering that in the t-AI diagram, the slopes of the 
incoming/outgoing radial light rays are each other's mir- 
ror images about vertical lines, it is evident from Fig. |21 
that AH± are spacelike everywhere except possibly in a 
neighbourhood of the central line M = and at future 
null infinity. From (|3.32|) and H3.30|l we have e = 8/3 and 
c = 2 in (|3.24|l . so we find i?AH- ^ 1 at the origin. Simi- 
larly by H3.31|l we have c — 3, and -Bah+ ~^ 1- Therefore 
AH^ are both spacelike at the origin too. 

Fig. |31also shows several outgoing radial null geodesies. 




FIG. 2: 3-d graph of the black hole formation process from 
Fig. Q the areal radius as a function of M and t. Each shell of 
constant mass evolves in a plane given hy M — const. It starts 
at i? = 0, then gets out of the past apparent horizon AH", 
then reaches maximum R, then falls into the future apparent 
horizon AH+, and finally hits the Big Crunch. Note that the 
surface intersects the R = plane perpendicularly all along 
the R — contour. The apparent horizons are intersections of 
the R{M, i) surface with the plane 7? = 2M. (Hence, it is easy 
to figure out how they would look in the k > Friedmann 
model, where ts = const, tc = const and all the ^(i)lM=const 
curves are indentical. 



not fulfilled, either the upper branch or the lower branch of some 
contours would be a non-monotonic function, whose derivative by 
M would change sign somewhere. At the changeover points, the 
tangents to the contours would be horizontal, and these would 
be the shell crossings. 



Each geodesic has a vertical tangent at the center. This 
is a consequence of using M as the radial coordinate. 
Since At/ AM = dR/dM/y/l + 2E on each geodesic and 
R oc M2/3 close to the center, so At/AM cx M~i/3 and 
At /AM — > oo as M 0. Each geodesic proceeds to 
higher values of R before it meets the apparent horizon 
AH+. At AH+, it is tangent to an i? = const contour, 
then proceeds toward smaller R values. The future event 
horizon consists of those radial null geodesies that ap- 
proach the AH+ asymptotically. In the figure, it lies be- 
tween the geodesies no 5 and 6, counted from the lower 
right corner of the figure; we shall discuss its location in 
more detail in sec. IIII F 31 

Geodesic no 5 from the lower right emanates from the 
center M = 0, where the Big Bang function has a lo- 
cal maximum. The tangent to the geodesic is horizontal 
there. This means that the observer receiving it sees the 
light infinitely redshifted, as in the Friedmann models. 
Geodesies to the right of this one all begin with a ver- 
tical tangent, which implies an infinite blueshift. These 
observations about redshift and blueshift were first made 
by Szekeres [s^. Likewise, the geodesies meet the Big 
Crunch with their tangents being vertical. 

By the time the crunch forms a.t t = tc(0), the future 
apparent horizon already exists (see Figs. nandOJ. The 
shells of progressively greater values of M first go through 
the AH, and then hit the singularity at t = tc{M). We 
assume that at the time t = t2, the singularity has al- 
ready accumulated the mass Ms, while the mass hid- 
den inside the apparent horizon at the same time is 
Mbh > Ms- Both of them grow with time, but at fixed 
t2 they are constants. From the definitions of Ms and 
Mbh it follows that 

t2 = tc{Ms) = tAH+{MBH)- (3.33) 

3. Location of the event horizon. 

Now we shall discuss the location of the event horizon 
in the spacetime model considered in sec. IIII Fl It will 
follow that, even though the model has a rather simple 
geometry, this is quite a complicated task that requires 
complete knowledge of the whole spacetime, including 
the null infinity. Hence, in a real Universe, where our 
knowledge is limited to a relatively small neighbourhood 
of our past light cone and our past worldline, and the 
knowledge is mostly incomplete and imprecise, the event 
horizon simply cannot be located by astronomical obser- 
vations. 

The future event horizon is formed by those null 
geodesies that fall into the future apparent horizon "as 
late as possible" , i.e. approach it asymptotically. Hence, 
in order to locate the event horizon, we must issue null 
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FIG. 3: Contours of constant i?-value (thinner lines) and outgoing radial null geodesies inscribed into the spacetime diagram 
of Fig. The _R-vaIues on consecutive contours differ by always the same amount. 
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geodesies backward in time from the "future end point" 
of the AH+. This cannot be done in the (M, coordi- 
nates used so far because the spacetime and the AH-I- 
are infinite. Hence, we must first compactify the space- 
time. The most convenient compactiflcation for consid- 
ering null geodesies at a null infinity is, theoretically, a 
Penrose transform because it spreads the null infinities 
into finite sets. However, in order to find a Penrose trans- 
formation, one must first choose null coordinates, and in 
the L-T model this has so far proven to be an impossible 
task, see Ref. 0, |^. Hence, we will use a less con- 
venient compactiflcation that will squeeze the null infini- 
ties into single points in the 2-dimensional (time - radius) 
spacetime diagram. It is provided by the transformation 



M — tan(/^), t = tan(T). 



(3.34) 



In these coordinates, the M} x space of Figs. ^"El 
becomes the finite [— 7r/2,7r/2] x [0, 7r/2] rectangle, see 
Fig. ^ The upper curve in the figure is the Big Crunch 
singularity; the future apparent horizon runs so close to 
it that it seems to coincide with it''. The horizontal line 
is the r = now line. The lower curve is the Big Bang and 
the past apparent horizon, again running one on top of 
the other. The point on the r-axis where the three lines 
meet is the image of the M ~ line of Fig. ^ squeezed 
here into a point because of the scale of this figure. 

The theoretical method to locate the future event hori- 
zon in Fig. 0] would now be to run a radial null geodesic 
backward in time from the point (//, r) = (7i'/2, 7r/2), i.e. 
from the image of the future end of the AH-I-. However, 
for the most part, the AH-I- runs so close to the crunch 
singularity, and the geodesies intersecting the AH-I- are so 
nearly tangent to AH+, that numerical instabilities crash 
any such geodesic into the singularity instantly, even if 
the initial point is chosen well away from (7r/2, 7r/2). This 
happens all the way down to /i = 0.5 at single precision 
and all the way down to /i = 1.1 at double precision. We 
did succeed, with double precision, only at ^ = 1.0, and a 
null geodesic could be traced from there to the center at 
/i = 0. At the scale of Fig. 01 this whole geodesic seems 
to coincide with the crunch and the AH-f . However, it 
is well visible if one closes in on the image of the area 
shown in Fig. ^ the closeup is shown in the inset. 

Actually, to make sure that we located the event hori- 
zon with an acceptable precision, we ran three differ- 



As can be verified from eq. 13.261 . the timo-difi^erence between 
the crunch and the AH-|- goes to infinity when M — > oo. How- 
ever, the ratio of this time-difference to the crunchtime goes to 
zero, which explains why the two curves in Fig. |l]meet at the 
image of the infinity. The same is true for the Big Bang and the 
AH-. 



ent null geodesies backward in time; one from the point 
/i = 1.0, r = To := rAi/+(1.0) right on the apparent hori- 
zon, another one from /i — 1.0, r = ri, where ti was in 
the middle between the AH-f and the crunch singularity, 
and the third one from /i = 1.0, r = r2, where T2 was be- 
low the AH+, with the time-difference tq — T2 = ti — tq. 
All three ran so close to each other that they actually 
coalesced along the way and reached the center as one 
curve; their coincidence in the inset in Fig. 0] is thus 
not an artefact of scale, but actual coincidence at double 
precision. 

Finally, the event horizon had to be transformed back 



1.5 



0.5 



-0.5 



-1.5 



AH+ & BC 




0.005 0.01 0.015 0.02 



AH- & BB 



0.2 



0.4 



1.2 



1.6 



FIG. 4: The spacetime diagram of Fig. 0compactified ac- 
cording to eq. (13.3411 . The whole region shown in Fig. 0is 
squeezed into the point where the three lines meet at the r- 
axis. The worldlines of the dust source would still be vertical 
straight lines here. The upper curve is the future apparent 
horizon and the Big Crunch singularity, they seem to coincide 
at the scale of this picture. The lower curve is the past ap- 
parent horizon and the Big Bang singularity, again coinciding 
only spuriously. The horizontal straight line is the r = now 
time. Inset; a closeup view of the image (in the coordinates 
(/.I, t)) of the region shown in Fig. The thicker line is the 
event horizon. It does not really hit the central point of the 
Big Bang; the apparent coincidence is just an artefact of the 
scale. More explanation in the text. 



12 



to the (M, t)- coordinates and written into the frame of 
the figuresn-|31 This is done in Fig. [S] As stated earher, 
the event horizon is located between the geodesies no 5 
and 6 from the right in Fig. O By accident (caused by 
our choice of numerical values in this example), the EH 
hits the center very close to the central point of the Big 
Bang, but does not coincide with it. 

This whole construction should make it evident that 
there is no chance to locate the event horizon by astro- 
nomical observations, even approximately. It only makes 
sense, in the observational context, to speak about an up- 
per limit on the mass inside the apparent horizon. This 
is why we identified the observed mass of a black hole 
with the AH in sec. IIVDP . 



5.14 




4.94 I 1 1 1 1 1 1 1 1 

0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 

FIG. 5: The event horizon (thicker line) written into the 
frame of Fig. Its intersection with the M = axis does 
not coincide with the central point of the Big Bang, this is 
only an illusion created by the scale. 



IV. THE GALAXY PLUS BLACK HOLE 
FORMATION MODEL 

Our aim is to model the formation of a galaxy with a 
central black hole, starting from an initial fluctuation at 
recombination. Our model consists of two parts joined 
together across a comoving boundary M = Mbh, with 
Mbh the estimated present day mass inside the black 
hole horizon. In the exterior part, we take existing ob- 
servational data for the present day density profile, and 
the initial fluctuation is made compatible with CMB ob- 
servations. For the interior, no observational constraints 
exist, so we propose a couple of possible descriptions, as 
detailed below. These are both LT models, and represent 
a collapsing body, and a dense Kruskal-Szekeres worm- 
hole in the sense of |27I. 



A. The black hole interior 

Astronomical observations do not say anything about 
that portion of galactic matter that had already fallen 
inside the apparent horizon by the time the electromag- 
netic signal that would reach the observer was emitted. 
What can be seen in the sky are only electromagnetic 
waves emitted by objects that were still outside the AH 
at the time of emission. Consequently, we are not con- 
strained in any way in choosing a model for the matter in 
the interior of a black hole, except for the need to match 
it smoothly to a galaxy model. 

The term 'black hole' is used in two disctinct ways. 
Firstly, there is the Schwarzschild-Kruskal-Szekeres black 
hole, which has the topology of two universes joined by a 
temporary wormhole, and begins its life as a white hole. 
This was generalised to a matter-filled version in [j^ . 
Secondly, there is the black hole formed by the collapse of 
a massive body, which has an ordinary topology without 
a wormhole. Only in their late stages (after the closure of 
the wormhole) do these two become essentially the same. 
Both of these can be reproduced by an L-T model, and 
we consider them in turn below. 

At this point, we must make a digression about the 
relation between model black holes such as those con- 
sidered here, and real collapsed objects that are called 
"black holes" by astronomers. There are two important 
points to be remembered: 



Note that this model in which the event horizon has been found, 
will actually only be used for the interior of our galaxy model 
— i.e. for that part of the galaxy that is invisible to outside ob- 
servers because of its proximity to the black hole horizon. The 
model of the visible part of the galaxy will be different, and to 



locate the event horizon there, the whole construction (integrat- 
ing a null geodesic backward in time from future null infinity) 
would have to be repeated. But, of course, the model of a single 
galaxy does not extend to infinity. Hence, it hardly makes sense 
to even speak of an event horizon in this context. 
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1. The matter proceeding toward a black hole disap- 
pears from the field of view of any real observer be- 
fore it hits the horizon, whatever horizon is meant 
(see below). For example, the event horizon is the 
boundary of the field of view for an observer at 
future null infinity, i.e. one who is infinitely dis- 
tant from the black hole and infinitely far into the 
future. Therefore, the observationally determined 
"mass of a black hole" is in fact only an upper limit 
of the mass that has actually fallen within the hori- 
zon; the latter can never be measured in reality. 
Models of black holes allow us to calculate better 
estimates of that mass, and even if the arithmetic 
difference between a model calculation and an ob- 
servational limit is small, it is important to under- 
stand the conceptual difference. 

2. ft is incorrect to speak about the event horizon in 
the context of observations. We may know where 
the event horizon is only in a model. In practice, 
we would have to be able to take into account the 
future fate of every piece of matter, including those 
pieces that have been outside our field of view up to 
now — an obviously impossible task. Worse still, 
if the real Universe is to recoUapse in the future, 
we might already be inside the event horizon and 
will never see any signature of it. Hence, the hori- 
zons whose signatures we have any chance to see 
(like disappearance of matter from sight, or a large 
mass being contained in a small volume) are appar- 
ent horizons — they are local entities, detectable 
in principle at any instant (although with the dif- 
ficulty mentioned in point 1). Even in our simple 
model considered in the previous section, it was 
rather difficult to determine the position of the 
event horizon, and it could be done only numeri- 
cally. 

In the following we will identify the estimated black 
hole mass, obtained by fitting a model to present 
day observations, with Mbh, the mass within the 
apparent horizon at time t2 (see eq. (|3.33|) ). 

B. A collapsed body 



We will assume that tc{M) is an increasing function, to 
be able to create a black hole, and to prevent shell cross- 
ings we assume t'g{M) < in addition. The functions 
have already been chosen so that the origin conditions 
are satisfied. 

To assure a smooth match to the exterior model, we 
require the continuity of the L-T arbitrary functions and 
their derivatives^ at M = Mbh- Given equations (|3.30|l 
and (|3.32|) . we solve for the constants a, b, T tso, 
so that E, ts, dE/dM & dts/dM are matched at the 
boundary: 



tso 
b 

To 



M dtB 
^^"TdM 



1 dtB 
2MdM 

M dtB 



M=Mbh 
M=Mbh 

4ttM 



2ttM^ dE 



6 dM 3(-2£;)3/2 (_2£;)5/2dAf 
1 dtB . 



27r 



3M2 dM 3M^{-2Ey/^ 



27r 



dE 



M{-2Ef/^ dM 



(4.1a) 
(4.1b) 



M=Mbh 

(4.1c) 



(4.1d) 



M=Me 



C. A wormhole 

Since we have no way of knowing anything about the 
matter and spacetime interior to AIbh, we can equally 
well fit in a dust-filled wormhole of the Kruskal-Szekeres 
type, constructed with the L-T metric. The essential re- 
quirement is that, at the middle of the wormhole, M must 
have a minimum value Mmin, and i?(A/min) = ~l/2. The 
minimum lifetime (time from past to future singularity) 
of the wormhole is then 27rA/min- We choose the following 
functions: 



tB = tBo - b{M - M^inf 



E = 



-M ■ 
2M 



+ a(M - M,ni„). 



(4.2a) 
(4.2b) 



For this example of the interior, we use the model of 
section UlI Fl but with different values of the parameters. 

® The Darmois junction conditions for matching an L-T model to 
itself across a comoving (constant r) surface only require the 



From these, the conditions for matching to an exterior at 
some given M value are 

matching of E and ts at the same M. 
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Mr, 



dAf I V \dM 1 



2dM 2M\ V \dM J 



M=Me 



M=Me 



dAf / I \ dM y \dM 



M dtB dE / ^ f dE\ 




J M=Mbh 



(4.3a) 
(4.3b) 
(4.3c) 
(4.3d) 




0.0017385272 - 



0.0017385270 



0.0017385268 



0.0017385266 



0.015 0.02 0.025 0.03 
M 




FIG. 6; The L-T model for the wormhole interior of section H V CI with parameters 15.211 . in geometric units. The second graph 
shows the neck of the wormhole magnified about 10^ times. Only on this scale are the apparent horizons distinguishable from 
the bang & crunch. The bang does actually curve downwards, but its variation is much less than that of the crunch. The neck 
is at the left where the mass reaches a minimum (Mmin = Mq). The apparent horizons cross in the neck at the moment of 
maximum expansion. The second sheet is not shown. It could be a mirror image of this sheet, or it could be quite different. 
The boundary with the exterior model (at Mbh ~ 0.03) is defined to be where t2 — tAH- Recombination ti is indistinguishable 
from the bang in the first graph, and far to the future in the second. 



This model was chosen for simplicity, and so is not very 
flexible. The matching fixes the value of A/min, which 
determines the lifetime of the wormhole. A model with 
more parameters would allow the wormhole lifetime to be 
a free parameter. The apparent horizons and singulari- 
ties for this model with the parameters (j5.2|l used below 
are illustrated in Fig. 



D. The exterior galaxy model 

1. The final density profile 

As our example for the density profile of the final state, 
we choose the galaxy M87. It is believed to contain a 
large black hole -i^l , and the density profile for its outer 
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part had been proposed some time ago [42 



=Po /(l + 5s' + cs* + ds6)" 



(4.4) 



where po = l.O-lQ-^s g/cm^, b = 0.9724, c = S.SIO-IQ-^, 
d = 2.753 • 10^^, n = 0.59. The distance from the center, 
s, is measured in arcmin, i.e. it is dimensionless, and is 
related to the actual distance r by 



r 21600 

D IT 



r6, 



(4.5) 



where D is the distance from the Sun to the galaxy. (Of 
course, no galaxy and no real black hole is spherically 
symmetric, so we cannot model any actual galaxy with 
the L-T solution. However, we wish to make our illus- 
trative example as close to reality as possible, and this is 
why we stick to an actual object.) For our purposes, we 
need density expressed as a function of mass, and a profile 
that goes to infinity at r — > 0, to allow for the singularity 
inside the black hole^*^. The mass profile corresponding 
to H4.4|) is not an elementary function. However, it turns 
out that the following very simple profile is a close ap- 
proximation to H4.5|l in the region considered in Ref. |43| 
(see Fig. [TJ: 



p{r) = po/('5r)'/', 



(4.6) 



with the same value of pq. The corresponding mass dis- 
tribution is 



M(r) = i^^por5/3<S-4/3^ 
5 



(4.7) 



and we make it singular at r = by adding a constant 
Ms to M{r) (this is the same Ms as defined before — 
the mass that had already fallen into the singularity at 
t2', so far this is still an arbitrary constant): 

M{r) = Ms + —npor'^^^S-^/^. (4.8) 
5 

Hence, r = [5(M - Ms)S'^^^ /12ttpo]^/'^ , and 

p{M) = (12^/5)4/Vo^/^'5~''^^(M - Ms)-^/^ (4.9) 

In principle, the value of p near to M — Mbh should be 
mcasureable, but in this regime the difference between 
the Newtonian and L-T definitions of density becomes 
too pronounced (because of the non-flat geometry). To 
infer p{Mbh) in a sensible way, the results of observa- 
tions should be consistently reinterpreted within the L-T 



Even though we will use the interior model near the centre, this 
requirement assists in joining interior and exterior smoothly. 



scheme, and such results are not, and will not be avail- 
able for a long time. Consequently, we will have to give 
up on this bit of information. 

From (|2.17|) . we find the corresponding R{M) via — 



(3/47r) ///^ Ax/p{x), which is 
/ ^ \3/5 

R{M) = — (5^/5 ^^^j „ ^^)3/5 ^^^^Q^ 

V127rpo/ 

Now we can determine the constant Ms by the require- 
ment that R = 2M at M = Mbh, i-c 



/ 5 



3/5 



Vl27rpo 
From here we find 



S^^' {Mbh - Ms) 



3/5 _ 



2M, 



BH- 



Ms = Mbh - ^ {2MBHf' 



5(54/3 



(4.11) 



(4.12) 



It follows, as it should, that Ms < Mbh- However, this 
result makes sense only if the Ms thus defined is positive. 
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FIG. 7: Comparison of the profile l|4.4^ (marked as r{x)) 
with 114.611 (marked as g{x)). The p on the vertical axis is 
mass density in units of 10~^^ g/cm'^, the x on the horizontal 
axis is the distance from the center in arc minutes, just as in 
Ref. [4^ . The values of parameters are given below eq. 14.41 . 
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The condition Ms > is equivalent to 

3/2 



24 



2M 



BH 



.52 



< 1. 



(4.13) 



For checking this inequahty, aU quantities have to be ex- 
pressed in geometric units. The black hole in M87 is 
beheved to have mass Mbh = 3 • lO^M© [Uliil, and 
its distance from the Sun is I? = 43 • 10^ light years. In 
geometric units, with 1 year = 31 557 600 s, c = 3 • 10^ 
cm/s, G = 6.6726-10-8cmVg-s2 andM© = 1.989-1033 g, 
this makes Mbh ^ 4.424- lO^^cm, S = 5.067- lO-^^cm"!, 
Pq = 0.741 - 10"53cin-2^ and the left-hand side of (jlT^ 
comes out to be 4.075 - lO^^i^ which is very safely within 
the limit. 



2. The initial fluctuation 

In order to define a model uniquely, we only need one 
more profile for the region M > Mbh, e.g. a density 
or velocity profile a,t t — ti < t2 or a specific choice 
of E{M) or tB{M). Apart from the density profile at 
t = t2, there are no other observational constraints in the 
region M > Mbh — the time by which galaxies started 
forming is not well known, and presumably different for 
each galaxy, nothing is known about the initial density 
or velocity distribution in the proto-galaxy at that time. 

Since the only quantities that are to some degree con- 
strained by the observations are the density and velocity 
profiles at the recombination epoch, it will be most nat- 
ural to use these for ti. Even so there is a problem: no 
numerical data are available for amplitudes of the tem- 
perature fluctuations of the CMB radiation at such small 
scales. We consequently chose a zero velocity fluctuation 
for one case and a zero density fluctuation for another. 
However it turned out that the latter was not suitable, 
since the solution required a collapsing hyperbolic region 
near Mbh in the exterior model. 



V. NUMERICAL EVOLUTION OF THE 
MODELS 

The programs written for papers I & II were adapted to 
facilitate this two-step model construction. First the ex- 
terior profiles were used as input to solve numerically for 
E{M) and tB{M) for Mbh < M < 1. The values of E k 
Ib and their derivatives at Mbh were extracted, and the 
parameters of the interior model calculated from them. 
Then the functions E h Ib were numerically extended 
into the interior model, down to M — or M — M^un- 
From this data, the model evolution was reconstructed 
using existing programs. 



Our first model uses the final density profile of section 
II V Dl for the galaxy at time t2 = 14 Gyr, and a flat ini- 
tial velocity profile at time ti = 10^ y, both exterior to 
Mbh- The interior of Mbh is a black hole formed by 
collapse, as described by H3.30|l and (|3.32|) . with parame- 
ters determined by the matching H4.1|) . Geometric units 
were chosen such that 10^^ Mq is the mass unit. In these 
units, the parameters are: 

a = 9.0662 x lO^'' , b = 0.012409 , 

tBo = 0.0017385 , To = 8.7279 x lO" (5.1) 

The resulting arbitrary functions and the behaviour of 
the combined model are shown in Figs |H1 and El Notice 
that the fluctuations of both density and velocity at re- 
combination are well within 3 x 10^'"' and 10^"*. The 
black hole singularity forms at time To = 13.618 Gyr 
(since iso = 2.7126 x 10"^ y is negligible), so it is 400 
million years old by today. 

Our second model uses the identical exterior, but the 
interior is a full Kruskal-Szekeres type black hole con- 
taining a temporary 'wormhole', as described by l|4.2|l . 
with parameters determined by the matching H4.3|l . The 
same geometric units were used, and, using the ' — ' sign 
in H4.3|l . the parameters are: 



a = -4.7475 x 10" 



b = 0.012409 



tso = 0.0017385 , M„ 



2.0562 X 10" 



(5.2) 



(Using the + sign in H4.2|l gives Mmin = 0.06 > Mbh = 
0.03, which is not acceptable.) Figs ^| and show 
the arbitrary functions and the behaviour of the com- 
bined model for this scenario. Because the exteriors 
are identical, the fluctuations of density and velocity at 
recombination (ti) outside Mbh are again well within 
CMB limits. The wormhole mass (minimum in M) is 
Mmin = 2.0562 Mq, and the future singularity first forms 
at Tq — 6.3613 X 10^^ sees after the past singularity. (The 
future and past black hole singularities are the extension 
of the crunch and bang into the middle of the wormhole.) 
The very short lifetime of the wormhole is a consequence 
of the need for E to go from —1/2 all the way up to 
— 1.7669 X 10~^ and arrive there with a negative gradient. 
(Even at constant M, igSl) impfies T{Mbh)/To > W^^, 
and the non-zero change in M only increases this factor.) 
Though models could no doubt be found with quite dif- 
ferent wormhole lifetimes, this example very effectively 
highlights the fact that the nature of the central black 
hole is essentially unknown. By recombination (ti), this 
black hole has accreted^^ 246200 Mq within the apparent 



The mass within the AH at ti is found by numerical root finding 
using 13.331 with ti instead of t2 and I3.26i 
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FIG. 8: The L-T model for the formation of a galaxy that develops a central black hole. Shown are the model-defining L-T 

functions E{M) & tsi^M), the pi[M) and 61 (M) fluctuations, the p2[M) and 62 (M) variations. The 62 (M) variation is zero, 
and only very small numerical error shows. Note that the graphs have been clipped at log(-R) = 0, log(6) = 15 and log(p) = —5 
(geometric units). The "EC" indicates the range considered is an elliptic region that is collapsing by ti. 



horizon, which is only 0.0048586 AU across. Any effect 
this might have on the CMB will not be observable for a 
long time. 



VI. SUMMARY AND CONCLUSIONS 



Wc have demonstrated the non-hncar evolution of 
an initial density perturbation at recombination into a 
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FIG. 9: The evolution of the L-T model for the formation of a galaxy that develops a central black hole. Shown are the 
evolution of R{t, M) and of p{t, M) for run ViOpfSO. In the R{t, M) graph, the origin R(t, 0) = is on the left, and expansion 
to a maximum occurs as time increases towards the "north-west". At the inital time, t\ = recombination, R> except at the 
origin. In the p{t, M) graph the view has been rotated by 180° relative to the R{t, M) graph for clarity, so that the orgin is 
on the right and time increases towards the "south-east" . By t2 ~ today, the innermost region has collapsed to a black hole of 
mass 3 x 10® Mq. Note that the graphs have been clipped at log(i?) = and log(p) = —5 (geometric units). 



galaxy with a central black hole at the present day, using 
the spherically symmetric L-T model. This is an applica- 
tion of the methods developed in papers I & II, in which 
an initial and a final state — each a density profile or 
a velocity profile — can be used to derive the arbitrary 
functions of an L-T model that evolves from one to the 
other. To correctly describe this process, a relativistic 
approach is necessary because Newtonian models are in- 
adequate for the description of black holes and their use 
inevitably leads to conceptual inconsistencies and con- 
tradictions. The L-T model is ideal for this purpose, as 
it has both Schwarzschild and Robertson- Walker limits, 
and a single model can describe a cosmology containing 
a black hole. 

For the final state at t2 = 14 x 10^ years w today, 
we chose the model of the mass distribution in the M87 
galaxy used in astronomical literature, eqs. (|4.4ll and 
(|4.5|) . More exactly, we approximated this mass distri- 
bution by a more elementary function whose values do 
not differ much in the range of interest, eq. H4.5|l . so 
that p can be calculated as an elementary explicit func- 
tion of the mass within a sphere of radius r. M87 was 
chosen since it is believed to contain a large black hole 
around its center, and several of its parameters have been 
measured or calculated. The inital fluctuation, at the re- 
combination epoch {ti = 10^ years), was chosen to be 
consistent with limits from the CMB, even though the 
smallest scales currently observable are much larger than 
those relevant to galaxy formation. We assumed zero ini- 
tial velocity perturbation, i.e. a Friedmannian velocity 
profile. This was sufficient for a unique numerical iden- 



tification of an L-T model that evolves the given initial 
state into the given final state. The resulting evolution 
of the L-T model, was found to be entirely reasonable — 
the implied initial density amplitude was well within the 
observationally allowed limit of 10~^, and the model was 
elliptic and already recollapsing by ^2 in the whole range 
of interest. Assuming the presence of a central black hole 
today, these profiles were taken to be valid for the exte- 
rior of the horizon, and a black hole model was smoothly 
joined on as the interior. 

A theoretical model of a black hole must necessar- 
ily include the accompanying entities: the final singu- 
larity, the apparent horizon and (whenever appropriate) 
the event horizon. We first discussed the general prop- 
erties of the apparent horizon. It must necessarily exist 
in every L-T model. The future apparent horizon AH+ 
exists in every collapsing model, the past apparent hori- 
zon AH" exists in every expanding L-T model; the ex- 
panding and recollapsing model has both AHs that can 
intersect each other only if there exists a neck or belly 
at which R' ^ M' = E' ^ t'g ^ 0, E = -1/2. In ev- 
ery case, the AH is given by the equation R = 2M that 
implicitly defines a function t — tAH{i')- We discussed 
whether the AH can be timelike, null or spacelike while 
going in or out. It turned out that only two cases are 
excluded: outgoing timelike for AH+ and ingoing time- 
like for AH~. The condition for the AH to be nontime- 
like is iAtAn/'^r < dM/dr, where € = ±1 for AH^ and 
A4 is the sum of rest masses of particles within the r- 
sphere. For the Fricdmann models, the AH is timelike 
everywhere. 
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FIG. 10: The L-T model for the formation of a galaxy around a pre-existing central black hole. Shown are the model-defining 
L-T functions E{M) & ts(M), the pi(M) and 61 (M) fluctuations, the P2{M) and 62 (Af) variations. (Again the 62 (M) variation 
is pure numerical error.) Note that the graphs have been clipped at log(J?) = 0, log(6) = 15 and log(p) = —5 (geometric units). 



The event horizon (EH) docs not always exist; for ex- 
ample it is absent in the A = Friedmann models. When 
it exists, the EH+ (EH~) is the collection of those future- 
directed outgoing (ingoing) radial null geodesies that 



approach the A11+ (AH^) asymptotically as i ^ +00 
(t —00). The equation of null geodesies is in general 
intractable by exact methods, and so locating the EH 
inevitably requires numerical integration. This can only 



20 




FIG. 11: The evolution of the L-T model for the formation of a galaxy around a pre-existing central black hole. Shown are 
the evolution of R(t, M) and of p{t, AI) for run ViOpfSO. In the R{t, AI) graph, on the left there is no origin, rather a wormholc 
with a tiny lifetime re-collapses to zero size. The thin flat wedge along the left is the growing singularity. (The apparent ripples 
are small numerical variations overemphasised by the graphics program.) In the p{t, M) graph, again with the view rotated by 
180° relative to the R(t, M) graph, the density diverges towards the singularity, and the flat wedge represents the part that has 
already collapsed. By t2 = today, the black hole mass has increased to 3 x 10^ Mq. Note that the graphs have been clipped 
at log(Ji) = and log(p) = —5 (geometric units). 



be (lone case by case, for specific forms of the L-T func- 
tions E{r) and tB{r). With a space of infinite volume, 
the numerical identification of the EH can most easily 
be done after the spacetime is compactified so that the 
future (past) edge of the AH+ (AH~ ) has finite time and 
radial coordinates. 

Then we calculated all the characteristic quantities of 
a black hole (the past and future singularity, both AH~s, 
radial null geodesies and the EH-I-) in a simple illustra- 
tive toy model with E <0 and displayed them in space- 
time diagrams. The model is recollapsing, has infinite 
total mass and volume, and a duration between bang 
and crunch that is finite at every finite mass, but goes to 
infinity as M — > oo and E ^ Q. 

Since no observational data exist (and, presumably, 
will not exist for a long time) concerning the interior 
of the horizon, two distinct forms of this central black 
hole — both L-T models — were considered, firstly a 
condensation that collapses to a singularity, and sec- 
ondly a full Schwarzschild-Kruskal-Szekeres type worm- 
hole topology. The parameters of these models were de- 
termined by matching their L-T functions to those of 
the exterior galaxy-forming model. In the case of col- 
lapse to a black hole, the central singularity is ~ 4 x 10^ 
years old by now. For the wormholc case, the final sin- 
gularity forms almost immediately after the Big Bang 
(6.36 X 10~^ sec), and is by today about as old as the 
matter in the galaxy. In this model, the black hole ac- 
cretes mass very fast into a very small volume, so that 
by recombination it had swallowed up 246 380 Mq in a 



region of diameter 0.00486 AU. However, all the numbers 
are strongly model-dependent, and there are no reliable 
observational constraints for model selection. The initial 
black hole is too small to have an observable effect on the 
CMB. Therefore we find both types of black hole are pos- 
sible. Perhaps small black holes, that avoid evaporation 
by rapid accretion, may seed galaxy formation. 

In fact, existing observational data does not have the 
resolution to constrain the initial data for our model. 
For example, as shown in Paper I, the perturbations of 
isotropy of the CMB temperature corresponding to single 
galaxies should have the angular size of « 4 x 10~^ de- 
grees, while the most precise current measurements have 
the resolution of 0.2°. 

Then, in the vicinity of the apparent horizon, 
the geometry of spacetime becomes measureably non- 
Minkowskian, while all the observational data available 
on mass distribution within galaxies were calculated by 
purely Euclidean reduction methods. We also stressed 
that what can be inferred from observations is only the 
upper limit of the mass inside the apparent horizon at a 
given time. It does not make sense to even speak of an 
event horizon in the observational context. 

In view of the paucity of data, our approach was the 
first exploratory step into an uncharted territory rather 
than an actual model to be compared with observations. 

The main limitation of the spherically symmetric L-T 
model is the lack of rotation, which slows collapse and 
stablises structures. Thus the model is good for much 
of the evolution into the non-linear regime, but becomes 
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less realistic as collapse sets in. 

Our results show that the L-T model is a very use- 
ful tool for this kind of investigation. However, for its 
parameters to be fine-tuned to results of observations, 
the observational data would have to be re-interpreted 
against the background of the L-T geometries. 
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Hence, in the Friedmann limit both branches of the AH 
are entirely timelike (outgoing in the expansion phase, 
incoming in the collapse phase) and monotonic with r. 
The derivatives and t'^ seem to become infinite at 

r = \j\fk. This is a coordinate effect. As seen from 
the metric ds^ = di^ - S'^it) (j^^ + r^dw^^ , there is 

a coordinate singularity at r = 1/^/k. Both derivatives 
become finite when the coordinates are changed so that 
e.g. r — sinr'. The quantity B in eq. (|A3|I does not 
depend on the choice of r. 

For a completely general Robertson- Walker model (i.e 
not just p = 0), repeating the whole reasoning, we obtain 
for the slope: 



APPENDIX A: THE FRIEDMANN LIMIT 

We briefly specialize the above results to the Fried- 
mann limit, where 

M = Mor^ 2E = -kr'^, 

is = const, R = rS{t), (Al) 

Mq > and fc > being arbitrary constants, and S{t) 
being the scale factor. The apparent horizon, where R — 
2M, has the equation 

S{t) = 2Mor'^. (A2) 
so by (Els, (ESIl and we obtain 

2M' 



B = l- 



R' 



S 



\/\ ~ kr'^ 



AH± 



(A3) 
(A4) 

(A5) 



SS 

which, after making use of the A = 
is equivalent to 



B 



-2p 



3p + P 



(A6) 

Einstein equations, 



(A7) 



This makes the AHs timelike for 1/3 > p/ p > —1/3, but 
spacelike for p/p > 1/3, or for p/p < —1. (Note that 
-6 = requires p = or divergent p.) Hence, it can- 
not be decided whether the AH is timelike or not with- 
out knowing the precise shape of the function S(t) (i.e. 
knowing the equation of state). It can only be said that 
as long as the source in the Einstein equations is ordinary 
matter known from laboratory (no cosmological constant 
or other self- accelerating medium), we will have S* < 0. 
Consequently, by (|A6|I . B < 0, which means that the AH 
will be outgoing in the expansion phase and incoming in 
the collapse phase. 
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